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A frame of a cone C is a minimal set of generators, and the lineality space L of C is the greatest 
linear subspace contained in C. Algorithms are described for determining a frame and the lineality 
space of a cone C{S) spanned by a finite set S. These algorithms can be used for determining the 
vertices, edges, and other faces of low dimension of the convex hull of a finite set H(S). All algo- 
rithms are based on the simplex method of linear programming. The problem of finding the lineality 
space can be successively reduced to problems in spaces of lower dimensions. 
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Let 



Introduction 



S = {A U . • ., A n } 



be a family of points in R'". We denote by L{S) its 
linear hull, by H(S) its convex hull, and by C(S) its 
conical or positive hull, i.e., the convex polyhedral 
cone expressible as the set of all nonnegative weighted 
sums of elements of S. We use the conventions 
L(</>)=C((/>)={0}. 

A subfamily T C S is called a frame [4J l of C(S) 
if C(T) = C(S) but C(T-{Aj})*C(T) for each AjE T. 
The greatest linear subspace contained in C(S) is 
called its lineality space [5]. Two main problems are, 
given S, to find a frame of C(S) and to determine the 
lineality space of C{S). These two problems are 
closely related. 

Several important problems are equivalent to or in- 
cluded in these two main problems. Consider, for 
instance, the system of linear inequalities 



AJX 



b h 



J- 



1, 



for XE: R'". Removing redundant constraints amounts 
to finding a frame for the cone spanned by the vectors 

Aj = l /M, j= 1, . . .. //, and I J. Determining the 
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dimension of a polyhedron given by linear inequalities, 
can be reduced to determining the dimension of a 
polyhedral cone C = {U\A r U ^0} (see [1()|). As to 
the latter problem, it suffices to note that the linear 
hull L{C) is the orthogonal complement of the lineality 
space of the polar cone C { ' = C{A\, . . ., A,,}, where 
A\, . . ., A n are the columns of A. This problem 
also arises if the uniqueness of an optimal solution to 
a linear program is to be established in the presence of 
degeneracy. 

Clearly, the problem of finding the vertices of a 
convex hull H(S) can be solved by finding the frame 
of a suitable cone. It should not be confused with the 
problem of finding the vertices of a polyhedron de- 
fined by linear inequalities. The latter problem cor- 
responds to finding the facets, that is, the proper faces 
of highest dimension, of some H(S). Even a moder- 
ate number of inequalities is apt to generate a huge 
number of vertices [6, 9], and by the same token, the 
number of facets of H{S) may be extremely large 
compared to the cardinality of S. Thus finding the 
facets of H{S) is inherently more difficult than finding 
the vertices of H(S). 

The problem of finding the edges of //(S), and other 
facets of low dimension, may still be expected to be 
essentially easier than to find the facets of H(S). In- 
deed, one is tempted to conjecture that an efficient 
determination of the facets requires prior determina- 
tion of the lower dimensional faces. 

Note that finding the edges of H(S) could be ac- 
complished by finding the vertices of H{S>2), where 

S2:={A m = UAi + Aj)\i<j}* 



Ri 



in brackel 



eferenees at the end of this paper. 



stands for "is defined by." 



However, this procedure is not recommended since 
the cardinality of S 2 tends to become quite large. We 
shall prefer a more compact technique, which works, 
with S rather than S-> or its analogs for higher dimen- 
sional faces, and which is closely related to an algo- 
rithm for solving the second main problem of deter- 
mining the lineality space of a cone. The faces of 
H(S) will be characterized in terms of sign patterns 
of matrices representing H(S). This aspect is ex- 
amined more closely in [13J. 

Charnes [2], Motzkin, Raiffa, Thompson, and Thrall 
[11], and Farrell and Fieldhouse [14] have proposed 
methods for finding the facets and the vertices of a 
convex polyhedron P determined by a system of linear 
inequalities. On the other hand, Goldstein [7] suggests 
an algorithm for solving the same problem where P is 
given as the convex hull of a finite set of points. 
Recently, Thompson, Tonge, and Zionts [12] published 
a method for removing redundant constraints. One 
of our algorithms for determining a frame of the cone 
C(S) essentially coincides with theirs. 

We shall give algorithms both for finding a frame 
for C(5), and for determining the lineality space of 
C(S). As will be pointed out, the latter algorithm 
can be modified to decide whether a given subset of 
S characterizes a face of C{S). Applied to single 
elements of S, such an algorithm can also be used to 
find the frame of a pointed cone C(S) and, therefore, 
the vertices of a convex hull H(S). 

All algorithms in this paper are based on the simplex 
method for linear programming (see, for instance, [3]). 
Not surprisingly, there are "primal" and "dual" 
algorithms in each case. We shall describe a primal 
algorithm for finding a frame and a dual algorithm for 
determining the lineality space of C(S). The remain- 
ing two algorithms are readily constructed by the 
reader. 

It has been our experience with the above algorithms 
that in the presence of degeneracies round-off errors 
may cause cycling if zero tolerances are not chosen 
and handled correctly. This danger is diminished if 
degeneracy provisions are made, but we shall not 
include such provisions in the descriptions of our 
algorithms. 



1. General Remarks 

A set S may contain more than one frame for C(S), 
and the cardinalities of these frames may differ. For 
instance, if C(S) is a linear subspace of dimension /, 
then the cardinality of a frame T may range from 
/+1 to 21. We shall only concern ourselves with 
finding some frame for C(S), and not, for instance, 
a frame of minimum cardinality. 

A point A c E S is redundant in S if C(S — {A C }) 
= C(S): otherwise, A c is necessary in S. The algo- 
rithms for finding a frame will consist of repeated 
applications of a subalgorithm which determines for 

some A c E S C S whether it is redundant or necessary 



in S. The subset S arises from S by deleting those col- 
umns that were already found to be redundant. 

Now let L be the lineality space of C(S). Clearly, 

L = C(S)nc(-S). 

We partition the set S into two classes, its lineal part 
Sl:=SC\L and its conical part Sc - = S — Sl. Then 



(1.1) 



L = C(S L ). 



PROOF. S/,CL, hence LDC(Sl). To prove the 
other inclusion, suppose IGL. Then X(EC(S) and 
XGC(-S). Hence X = ZAni u m ^ 0, andX=2(-Advu 
vi 2* 0. It follows that 

If m + Vi > 0, and in particular if u\ > 0, then clearly 
-At E C(S) and, therefore, A { E S L . Thus X E C(S L )-' S . 
Let 5 be a subset of S such that S'DSi,. Let L 
denote the lineality space of C(S). Then (1.1) gives 

l = C(S0 n C(-Sj) c c(S) n q-s) c qs) n C(-S) = l , 



whence L = L. As a consequence 

(1.2) S L = S L for S L CSCS. 

We shall describe an algorithm for finding the lineal 
part Sl of a given set S. Again, this algorithm will 
be based on a subalgorithm, which decides whether 
a specified A c E S C S is in Sl- The set S arises 
from S by deleting those points A\ which were already 
found to be in Sc. This deletion is justified by (1.2). 

We proceed to show how the problem of finding 
faces of C(S) relates to the problem of determining 
lineality spaces. The faces of a cone C are extreme 
subsets, that is, if the sum of points in C lies on a face 
F of C, then so do all the summands. More precisely, 

(1.3) X E F \ and X = HYiiii with F, E C and u\ ^ 0, 
imply m = for Yi $ F. 

This property characterizes faces. If F is a face of 
C(5), and if S F ;=FHS, then 

(1.4) L(Sf) is the lineality space of C(S U (— Sf)). 

Proof. C(S U (-Sf)) C C(S f U (-S f )) = L(S f ). 
Now suppose 

x $ qs u (-s F )) n C((-S) u s F ). 

Then there exist U[, wu vu Zi ^ such that 

X= ^ Am- 2 Am=- 2 A > v >^ 2 A ^- 

A\ E S A, E S F A] E S A; E Sr 



3 The symbol "— " marks the end of a. proof. 



Hence 



2. Cones and Equivalent Matrices 



2 Aj(m + V})= ^T Ai{iv,-\- zi). By definition, the point A c is redundant in S if and 

A,es /e.sv only if C{S — {A c }) = C(S). This in turn is equivalent 

to A c G C(S-{A C }), or 

By (1.3), iij-\-Vi = if Aj^Sy. In particular, Ui = 

if /*, ^ S,. Hence (2 .1) ,4 C = 2 4 a * ^' * °- 



X = 2 4m- ^ ^ 



which implies /¥ E L{Sf). — 

The points 4, . . ., A k (all in S) are said to s*/6- 
determine a face F of C(S) if ^ ' 



The point A c is in the lineal part Sl of 5, if and only 
if A c E L = C(S) H C(—S). This is equivalent to 
-^ c G C(S-{/4 c }), or 



■A c ="2AjV h vj^O. 



F = L{A U . . . ,4-} flC(S). 



We then have 

(1.5) 77ie points Ai, . . . , A k i/i S sab determine a face 
of C(S) i/ and on/y i/" L{Ai, . . . , A k } is the 

lineality space of C(S U {— Aj, . . . , — Aa}). 

Proof. Suppose L{A U . . ., A k } is the lineality 
space of C(SU {—Au • • •. — ^a}). In order to 
prove that F:=L{Au • • ., A k } D C(S) is a face, 
we have to show according to (1.3) that if X E F, 
X=XAjUj with Uj ^0, and ///, >0, then An E F. Sup- 
pose, therefore, that X and iin are as above. Then 



A h =— (Z + Z), 



/here 



Zf-^^GQ(-.S)UR, . . ., 4t}). 



Now X is in L{A\ % . . ., /4a-}, which is by hypothesis 
the lineality space of C{S U {— A u ■ • • , —/4a-}). 
Therefore, X E C((-S) U {A u . . ., A}), whence 
A h E C((-S) U {A u . . ., ^a}). It follows that A h 
is in the lineality space of C(S U {— A u . . ., —A k }). 
Thus A h E L{A U . • ., 4} n C(S) = F. 

On the other hand, if F:=L{A U . . ., A k } f) C{S) 
is a face, then L(Sf) = L{A U . . ., A k } is the lineality 
space of C(S U (-S F )) = C(S U {-A u . . ., -/4a-}) 
by virtue of (1.4). — 

Note that if 4, . . ., A k subdetermine F, then 
L(F) = L{A U . . ., A k }. Hence 

dim (F) = A 

if A u • • -, /4a are linearly independent. Thus any 
algorithm for determining lineal parts can be used for 
deciding whether a given set of points {A u . . . , A k } 
C S subdetermines a A- -face (= face of dimension k) 
of C(S). 



The standard "phase /" simplex procedures of linear 
programming are available for deciding the solvability 
of (2.1) or (2.2) (see, for instance, [3]). All that re- 
mains is to integrate such a procedure with the overall 
algorithm in a way that minimizes the rearrangement 
of data. 

The points A u . . ., A n of S can be represented by 
an m X n matrix 



A:=(Ai 



A fl ). 



The matrix A is in canonical form if it contains an 
m X rn permutation matrix, called basis. Its columns 
are the basic columns; all other columns are nonbasic. 
If A is in canonical form, and if the column A; is non- 
basic, then the elements of Aj constitute the represen- 
tation of the point Aj in terms of the basis of A. Cri- 
teria for redundancy and lineality will be therefore 
particularly simple if thejnatrix A is in canonical form. 
Two matrices A and A are called equivalent if there 
exist (possibly rectangular) matrices 7" and T such that 
TA=A and TA=A. If and only if A and A are equiva- 
lent in this sense, then for all X E R": 

(2.3) AX = if and only if AX = 0. 

Both (2.1) and (2.2) express the existence of a suitable 
linear dependence 

AX = 0< X^O 

of the columns of A. By (2.3) equivalent matrices 
admit the same linear dependences. If, therefore, 
A is replaced by any equivalent matrix A— preferably 
in canonical form — then the same columns (= columns 
with the same indices) will be redundant in the re- 
spective column sets, and the same columns will be 
in the lineal parts. 



3. Determining a Frame of C{S) 



For matrices A in canonical 
criteria are immediate: 



form the following 



(3.1) If A } is nonbasic andAj^O, then A j is redundant. 

(3.2) If Aj is nonbasic and has exactly one positive 

entry a rj , then the basic column A c for which 
a rc = 1 is redundant. 

The following two criteria hold whether A is canonical 
or not: 

(3.3) If the \th row of A, henceforth denoted jA, con- 

tains exactly one negative entry a is , then A s 
ls necessary. 

(3.4) //" jA contains exactly one positive entry a is , 

£/ica? A s is necessary. 

As was pointed out before, any "phase /" simplex 
procedure can be used for deciding whether a given 
column is redundant. We prefer a variant (described 
in [8J) without artificial variables which works on each 
infeasible row separately, treating it to some extent 
as an objective function while conserving the feas- 
ibilities already achieved. While this variant may not 
be the most efficient one for finding a first feasible 
solution — using a positive combination of the infeasible 
rows is in general better — the terminal situations of 
this variant are precisely the ones to which criteria 
(3.1) and (3.3) apply; it is therefore particularly easy 
to implement. 

The matrix A will be repeatedly transformed and 
some of its columns may be deleted. In the following 
description of the algorithm, the symbol "v4" will 
always refer to the particular matrix at hand. It will 
also be convenient, not to change the indexing of 
columns which remain after others have been deleted. 



(3.5) 



Primal algorithm for determining a frame: 
(i) {Canonical form) Use Jordan elimination 
for bringing A into canonical form. 
Every column is labeled "undecided." 
(ii) (Constant column) Select a nonbasic 
undecided column A c as "constant col- 
lumn." If there are no such columns, 
go to (viii); else proceed to (hi), 
(iii) {Pilot row) If A c ^ 0, delete A c and re- 
turn to (ii). Else select pA such that 
a pc < 0; call it the "pilot row." 
(iv) (Pivot column) If a pc is the only negative 
entry in the pilot row, change the label 
of A c from "undecided" to "necessary" 
and go to (ii). Else select a "pivot 
column A i such that a p \ < and / t^ c. 
(v) {Pivot row) Select a "pivot row" r A 
such that 



= 



Ore 

a,i 



an 



0, an > 



There exists always an index r of this kind, since either 
the minimum on the right-hand side is finite and as- 
sumed—we use the convention min (</>) = + oc — - 5 r 
r=p is a permissible choice. 



(vi) (Pivoting) Execute a simplex step (= Jordan 
transformation) with a r i as pivot. In the 
absence of degeneracies, this will increase 
the old entry a pc while keeping nonnegative 
entries of A c nonnegative. 

(vii) (Return) If the new entry a pc is still negative, 
keep the pth row as pilot row and go to (iv). 
Else go to (iii). 

(viii) (Termination) If all basic columns are de- 
cided, terminate the procedure. Else select 
an undecided basic column A c . 
(ix) (Clear basis) Suppose a rc =l- If this is the 
only positive entry in r A, then change the 
label of the column A c from "undecided" to 
"necessary" and go to (viii). Else pivot so 
as to remove the undecided column A c from 
the basis. Go to (iii) with A c as constant 
column. (End of the algorithm.) 

In the absence of degeneracies, this algorithm in- 
creases at every step the entry a pc until it becomes 
nonnegative or remains the only negative entry in the 
pilot row. In the latter case, the column A c has been 
decided. In the former case, a pc stays nonnegative 
during all subsequent transformations, while some 
other negative element of A c is being increased. 
Hence a decision on Ac will be reached after finitely 
many simplex steps. 

The algorithm is sped up by checking after each 
iteration whether some of the columns or rows satisfy 
criteria (3.1), (3.2), (3.3), or (3.4). 

A dual algorithm results if the subalgorithm which 
decides whether A c is redundant or necessary is re- 
placed by its dual. The primal decision algorithm 
tries to make a given column of A nonnegative (cri- 
terion (3.1)), and encounters criterion (3.3) if this is 
not possible. The dual decision algorithm aims at 
criterion (3.3), that is, it attempts to make all entries 
but one in some row of A nonnegative. If this attempt 
fails, then criterion (3.2) is encountered. This algo- 
rithm is a "phase /" procedure for the dual simplex 
method. 



4. Determining the Lineality Space of C(S) 

The following criterion is immediate for matrices 
in canonical form: 



(4.1) 



//Aj ^0, then Aj E S, 



For any matrix A one has 

(4.2) //jA 5* 0, then all columns Aj with ay > belong 
to S c . 

An algorithm analogous to the algorithm (3.5) can 
be based on these criteria. The pivot rules arise 
from those of the algorithm (3.5) by reversing the sign 
of the constant column. Whenever a nonnegative 
row \A is found, all columns Aj with ay > are deleted. 



Thus, ith row of the remaining matrix vanishes and 
is therefore deleted. Whenever a nonpositive column 
A c is found, it is labeled "lineal," and this process is 
repeated until an empty matrix results or until all 
columns are labeled "lineal." There are primal and 
dual strategies available, the former aiming at non- 
positive columns, the latter at nonnegative rows. 

The algorithms are sped up if the following refine- 
ment of criterion (4.1) is used (for matrices in canonical 
form). (See [13]): 

(4.3) Suppose A c ^ 0, and let A be the matrix that 
arises from A ifjdl rows 4 A with a ic < are 
deleted. A Then Aj is in the lineal {conical) 
part of A if and only if Aj is in the lineal 
(conical) part of A. If A c < 0, then S c = 0. 

PROOF. Not only A c but also all basic columns Ah 
such that Oih = 1 for a\ c < are in S L . These basic 
columns Ah span the linear subspace E := {X\xi = if 
a ? > = 0}, which is therefore contained in the lineality 
space L of C{S). We define Aj by putting 




Then Aj — Aj E E C L. Hence Aj E L if and only if 
Aj EL. 

Denote by S the set of all Aj. It clearly suffices to 
show that Aj E Sl if and only if Aj E Sl- „ Suppose 
Aj E Sl. Then Aj E L and therefore — Aj = LAiUi, 
Ui ^ 0. This relation remains true^if we replace A\ 
by Af. —Aj = XAiUu whence Aj E Sl. In the other 
direction, if — Aj = ^AjUi, w^O, then 



- Aj = Aj — Aj + X(A j — A i)m + 2/1 tui 



L C C(S). Hence 



-/^•ECOS) 



but Aj — Aj, Aj — A 

and, therefore, /fj E S/,.— 

We proceed to describe a dual algorithm. The cor- 
responding primal algorithm is readily constructed by 
the reader. Again it will be convenient to denote by 
A the particular matrix at hand, and not to change the 
indices of rows and columns if other rows and columns 
are deleted. 

(4.4) Dual algorithm for determining the lineality 
space of 'C(S): 

(i) (Canonical form) Use Jordan elimination for 
bringing A into canonical form. Label all 
zero columns of A "lineal." 

(ii) (Termination) If A is empty or all columns 
of A are labeled "lineal," then terminate the 
procedure. Else proceed to (iii). 

(iii) {Objective row) Select any row C A as "objec- 
tive row." 

(iv) (Pilot column) If r A^Q, then delete all 
columns Aj with a C j > 0, delete V A , and return 
to (ii). Else select a "pilot column" A p 
with a C p < 0. 



(v) (Pivot row) If A p ^ 0, delete all rows -,A with 
a\ v < 0, label all zero columns generated by 
this deletion (all columns, if A,, < 0) "lineal," 
and return to (ii). Else select a "pivot row" 
,-A with a rp > and r ^ c. 

(vi) (Pivot column) Select a "pivot column" Ai 
such that 



a C j ^ 0, Orj < 



— 5* max \ — 

Or! [CLrj 



(compare algorithm (3.5) step (v)). 
(vii) (Pivoting) Pivot on a,i. In the absence of 
degeneracies, this will increase the entry 
a cp while keeping nonnegative entries of C A 
nonnegative. 
(viii) (Return) If the new entry a C p is still negative, 
keep the /;th column as pilot column and go 
to (v). Else return to (iv). (End of 
algorithm) 

In the absence of degeneracies, this algorithm will 
terminate in a finite number of steps (compare (3.5)). 
The algorithm is sped up by checking after each itera- 
tion whether some of the rows or columns satisfy cri- 
teria (4.1) or (4.2). 

5. Determining h -Faces of C{S) 

For finding all A -faces of C(S) we employ a subalgo- 
rithm which decides for each set of ATinearly independ- 
ent elements of S whether or not they subdetermine a 
face. We do not raise the question of the most effi- 
cient arrangement of this subalgorithm within the 
overall search algorithm. Instead, we try to formu- 
late the decision algorithm in a manner that does not 
preclude its implementation as a subalgorithm. To 
be more precise, in order to decide whether A u . • >,Au 
subdetermine a A*-face, it suffices by (1.5) to determine 
the lineality space of C(SU{— A u • ■ -, ~A k }). If, 
however, algorithm (4.4) is used for this purpose, 
after the columns — A\, . . .,—/*/> are adjoined to 
the matrix A, then a matrix results, which is no longer 
equivalent to A, and each decision must thus start 
from scratch. This procedure is therefore not suitable 
for implementation as a subalgorithm. 

We proceed to sketch an algorithm which does not 
require new columns to be adjoined, and which ter- 
minates with a matrix that is equivalent to the original 
one. 

A column Aj is lexico-nonnegative if <7ij>0, or 
u\j = and a 2 j > 0, or a\j = ci2j = and a$j > 0, or . . ., 
a\j= - • - = a„j = 0. If Aj is lexico-nonnegative and 
does not vanish, then Aj is lexico-positive. We call 
a matrix lexico-nonnegative (lexico-positive) if this 
holds for each column. We then have the following 
generalization of criterion (4.2): 

(5.1) Suppose A arises from A by deleting some rows. 
If A is lexico-nonnegative and Aj # 0, then 
Aj £ So 



PROOF. Let B be the row that results from com- 
bining the rows of A with weights e, e 2 , . . . . For 
suitable e > we have B ^ with bj > if and only if 

Aj #0. Then criterion (4.2) applies to the matrix ( R l, 

which is equivalent to A. 

We proceed to formulate two criteria on which an 
algorithm for finding the Maces of C{S) can be based. 
Suppose A u . . . , Ak are basic columns. We call 
the submatrix A of A formed by all those rows in which 
A u • • ., Ak simultaneously vanish the complement 
of A\, . . . , A/ x -. The complement of the entire basis 
of A is the empty matrix. Schematically, we have 



/(A") 










fim-k) 





A = 



Criterion (5.1) then leads to (see [13]): 



j=compl 



ement 



(5.2) 



// the 

A„ . 
A,, . 



complement of the basic columns 
., A k is lexico-nonnegative, then 
., A k subdetermine a k-face. 



Proof. It follows from (5.1) that A u . . ., A k , 
— Au . . ., —Ajc and those elements of S which are 
linear combinations of the former are the only elements 
in the lineal part of S U {—Au • . ., — A k }, and (1.5) 
applies. — 

(5.3) The basic columns Ai, . . ., A k do not subde- 
termine a face if their complement contains 
a nonpositive nonzero column. 

PROOF. Assume for simplicity that Ak+i, . . ., A,„ 
are the remaining basic columns, and that au=l for 
i^m. The complement of Au • • •-, A k then con- 
sists of the rows \A with i > k. Let now A c , c > m, be 
a column whose portion in the complement is nonposi- 
tive and nonzero. Then 



A c =^ Atut+ Y A\wu 



where Wi = a; c ^0, and, therefore, A c e C(( — S) U 
{A\, . . . , Ak}). Hence A c belongs to the linearity 
space of C(SU {-A u . . . , -A k }). Now if 
A u . . . , Ak subdetermine a face, then L{A U . . . , 
Ah} is the linearity space of C(S U { — A u . . . , 
-A k }) by (1.5). Thus A c e L{A U . . . , A k }, which 
contradicts the fact that not all w\ vanish in the above 
unique representation of A c in terms of the basis.— 
It is easy to see that algorithm (4.4) and its corre- 
sponding primal algorithm can be modified to yield 
an algorithm whose terminal situations are (5.2) and 
(5.3). No columns and no nonnegative rows are 
deleted. A hierarchy of nonegative rows leads even- 
tually to the nonnegative submatrix required by (5.2). 



The pivot columns vanish in these rows. Hence 
pivoting does not spoil the nonnegativities already 
achieved. 



6. The 1 -Skeleton of H(S) 

The problem of finding the vertices of a convex 
hull H{S) can be reduced to finding the frame of the 
coneC(5) where S \ = {A i = [A i \\A i eS}. Similarly, the 

edges of H(S) correspond to the 2-faces of C(S). There 
are, however, special features of the problem of finding 

the vertices and edges of H(S). First, 

(6.1) If A is equivalent to A, then Aj€C(S— {Aj}) 
if and only ifA s e H(S — { Aj}). 



Proof. If Aj = ^ Am, then by (2.3) 4/ = ]£ Am, 

which gives Zut = l by considering the last component. 
As a consequence of (6.1), a criterion for necessity 
can be based on the relative magnitudes rather than 
the signs of the elements of A (Goldstein [7]): 

(6.2) If A is equivalent to A, and if a^ is the unique 
maximum or minimum of the entries in a row 

jA of A, then Aj is a vertex of H(S). 

Indeed, Aj is a vertex of H(S), and by (6.1) necessary 
in S. Criterion (6.2) is valuable in practice since it 
enables one in general to find some vertices right away 
as well as during the algorithm. It can be generalized 
to the case in which the maximum or minimum is not 
unique. Suppose the maximum of the first row is 
assumed for the subset M C S. Then a column 
Aj e M is necessary in S if it is necessary in M. Hence 
the examination can be temporarily restricted to M. 
In particular, criterion (6.2) can be applied to M. 
For instance, if ao s is the unique maximum or minimum 
of the entries a 2 j with A } eM, then A s is necessary 
m S. 

Suppose that A is in canonical form and equivalent 
to A. Then the same argument that established (6.2) 
gives 

m 
i= 1 

for all / Thus each column has at least one positive 
entry. This follows also from the fact that C(S) is 
pointed i.e., has linearity dimension = (4.2) and that 
S contains no zero columns. 

Suppose S is a frame. Then each nonbasic column 
has at least one negative and two positive entries by 
(3.1) and (3.2). If a column has precisely two positive 
entries, then by (5.3) the corresponding basic columns 
do not determine an edge. 
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